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ABSTRACT 


Spatial Modal Testing is a parameter estimation approach, which provides spatial 
information regarding the dynamics of a structure, at a large number of points. The 
classical approach of obtaining frequency response function (FRF) estimates offers 
relatively little spatial information, for the emphasis is on observing the 
machine/component dynamics at a large number of frequencies at a limited number of 
specified locations on the machine/component. The FRF approach provides good 
accuracy on estimates of natural frequencies and damping ratios. The correctness of the 
estimates of mode shapes is, however, not comparable. Advances in sensor technologies 
facilitate measurements at a relatively large number of structure locations, and spatial 
modal analysis procedures seek to remove the deficiencies in the existing mode shape 
estimation procedures. The present study is an extension of some existing works on 
spatial modal analysis, dealing with single and two mode approximations. An algorithm, 
for estimation of the fundamental and desired number of successive modes, from spatial 
response data at multiple excitation frequencies, has been developed in the present study. 
Approximations of beam deflections, under specified harmonic excitation, are 
approximated through Forsythe polynomials, which can be compared with actual 
experimental measurements of the deflected beam. A Least Square technique is used to 
minimise the error between measured and Forsythe approximation data. Orthogonality 
properties of the Forsythe polynomials and normal mode functions are further employed 
to estimate individual mode shapes. Natural frequencies are obtained through standard 
relationships. An algorithm has been developed, for carrying out the procedure, 
successively, at various excitation frequencies, estimating two adjacent modes each time 
and simultaneously, refining existing previous modal information. The algorithm is 
demonstrated through computer simulations of the response of beams, with some 
standard boundary configurations. The analysis has been carried out in terms of non- 
dimensional parameters. The accuracy of the results is illustrated through comparisons of 
the estimated mode shapes with the ‘exact’ ones, available from classical beam theory. 
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CHAPTER 1 


INTRODUCTION 


Obtaining mathematical descriptions of the dynamic behaviour of machine 
components, is a major activity in vibration analysis. Modal Testing encompasses 
such experimental studies, for mathematical modeling. The form of the mathematical 
description or model varies considerably from one application to the other. It can be 
an estimate of natural frequency and damping factor, in one case, and estimating full 
mass-spring-dashpot model for another. One of the earliest developments in the area 
was by Kennedy and Pancu (1947) who evolved descriptions like ‘Resonance 
Testing’ and ‘Mechanical Impedance’, for applications involving natural frequency 
and damping determination of aircraft structures. Bishop and Gladwell (1963) 
provided methods for more precise measurements and applications. Major advances 
in transducers, electronics and digital analyzers, by 1970, helped to firmly establish 
modal testing as a powerful tool of experimental analysis. Reference can be made to 
papers by Ewins and Griffin (1981), Allemang (1984) and Mitchell (1984), for 
bibliographies of studies pertaining to basic modal testing methods. The text, by 
Ewins (1986) is widely referred to as a treatise encompassing most of the aspects of 
classical modal testing. 

Estimation of natural frequencies and description of mode shapes are the most 
common modal testing applications. These permit their model identification and allow 
correlation with a theoretical model. In a classical approach this is done through 
obtaining frequency response function (FRF) estimates from time dependent forces 
and responses spatially distributed about the structure. The denominator, of the 
frequency response function, comprises of the spatially invariant natural frequency 
and damping parameters, while the numerator contains the spatial information on the 
mode shape of the structure. 

The above classical procedure offers, relatively, little spatial information. For 

example, a typical test may have ten measurement locations on a structure, which is 
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provided harmonic excitation at a specified location. The excitation frequency is 
swept over the frequency range of interest and a Frequency Response Function (FRF) 
plot is obtained. These plots are then utilised for estimation of system parameters, like 
natural frequency, damping ratio, mode shape etc. While good accuracy is obtained 
on estimates of natural frequencies and damping ratios, the correctness of the 
estimates of mode shapes is not comparable. Subsequent generation of mass and 
stiffness matrix models is also not accurate. The primary reason, for these 
inaccuracies, is that response measurements are Hot made at sufficient number of 
locations on the structure. Advances in sensor technology and the advent of 
instruments like the Laser-Doppler Vibrometers, have facilitated measurements at a 
relatively large number of structure locations, in contrast to conventional sensors like 
accelerometers or proximity probes. Moreover, the measurement points can located 
with a much greater degree of accuracy, through a Laser beam measurement 
technique.The alternative to this approach is to obtain spatial information at many 
more points at a fixed frequency. This has recently become realizable through 
developments in Laser Doppler data acquisition techniques. The scanning laser is able 
to acquire velocity information at a number of spatial locations on the structure with 
sufficient density, accuracy and speed to offer the possibility of extraction of modal 
parameters from spatially intensive arrays as opposed to the classical frequency 
intensive format. The need for exploring the spatial data intensive procedure arises 
from the problems faced in validating analytical structural models. These models, 
typically finite element models, are characterized by a large number of spatial degrees 
of freedom. 

Reference can be made to the works of Archibald and Wicks (1991), Archibald 
(1992), and by West et al (1993), for attempts towards spatial modal analysis. 
Standard autoregressive models as in Pisarenko Harmonic Decomposition (Pisarenko, 
1973) and Oversize Total Least Squares, have been used to approximate a normal 
mode function of the beam and are shown by yield unbiased results in the presence of 
measurement noise. West et al (1993) have also, alternatively employed the recursive 
form of Forsythe polynomials to approximate a normal mode function. The 
orthogonality property of Fosythe polynomials, yields easier least square formulation 
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and removes the ill conditioning problem while extracting the Forsythe coefficients 
from measured data. 

The above attempts were, however, restricted to approximation of a single normal 
mode function of a beam by a polynomial function. The steady state vibration 
amplitude was obtained at various points along the length of a beam, for an excitation 
frequency, very close to one of its natural frequencies. The excitation frequency being 
very close to a natural frequency, it was assumed that contribution from adjacent 
modes is negligible and the deformed beam shape resembles the corresponding 
normal function of the beam. This is, however, not be always possible for, firstly it 
requires an a-priori knowledge of the natural frequencies of the beam and secondly, it 
may not always be possible to place the excitation force accurately enough on the 
structure, so as to constrain the response to a particular mode. 

Vaze (1998) developed a general procedure, for a beam excited at a specified 
frequency and location, to estimate the modal participation of adjacent modes, from 
steady state spatial data. He employed recursive form of Forsythe polynomials for his 
analysis. Characteristic properties of beam modes lying on either side of the excitation 
frequency were determined. The analysis was thus, restricted to estimation of two 
natural modes of the beam. 

The present study attempts to extend the work, cited above, for estimation of more 
than two natural modes. An algorithm, for estimation of the fundamental and desired 
number of successive modes, from spatial response data at multiple excitation 
frequencies, has been developed. The mode shapes of the beam have been 
approximated through Forsythe polynomials and an approximation of the beam 
deflection under a given excitation is generated. This approximation can be compared 
with actual experimental measurements of the deflected beam. A Least Square 
technique is used to minimise the error between measured and Forsythe 
approximation data, through adjustments of the coefficients of Forsythe polynomials. 
Orthogonality properties of the normal mode functions are further employed to 
decompose the Forsythe polynomial representation into the individual mode shapes of 
the beam contributing to the response. Natural frequencies can be obtained through 
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standard relationships. An algorithm has been developed, for carrying out the 
procedure successively, at various excitation frequencies, estimating two adjacent 
modes each time and simultaneously, refining existing previous modal information. 

The algorithm is demonstrated through computer simulations of the response of 
beams with some standard boundary configurations. Analysis has been carried out in 
terms of nondimensional parameters. Accuracy of the results is illustrated through 
comparisons of the estimated mode shapes with the ‘exact’ ones, available from 
classical beam theory. 

Chapter 2 describes the spatial modal analysis procedure. The algorithm for 
estimation of successive mode shapes along with the numerical studies is given in 
Chapter 3. Chapter 4 gives the conclusions of the present study and outlines the scope 
for future work. 
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CHAPTER 2 


ANALYSIS OF SPATIAL DATA 


The procedure, developed for analysis of steady state, forced vibration response data 
for estimation of modal participation, is discussed in this chapter. The beam is excited 
at a desired frequency and contribution, of the two natural modes of the beam, lying 
on either side of the excitation frequency, to beam response is estimated. The focus is 
on obtaining accurate descriptions of the mode shapes. Estimates of natural 
frequencies result as a by-product of the procedure. 

The beam response is obtained at a number of spatial points along its length, in terms 
of the in-phase and out-of-phase vibration amplitudes. Recursive form of Forsythe 
polynomials is employed to approximate the deflected beam shape. Least square 
technique is employed to minimise the error between the approximated and measured 
responses, to obtain the coefficients of Forsythe polynomials. Orthogonality 
properties of normal mode functions are further employed to decompose the Forsythe 
polynomial representation into individual mode shapes of the beam contributing to the 
response. Standard relations give the natural frequencies from estimated spatial 
frequencies. 

2.1 Forsythe Polynomials 

Orthogonal polynomials are most suited to approximate the dynamic deflection of 
beams, in view of the orthogonality property of their natural modes. Forsythe 
polynomials have been used by West et al(1993) to approximate the fundamental 
beam mode. They found that such a formulation of beam response is easier and 
removes ill conditioning. In view of these properties, Forsythe polynomials have been 
employed in the present study also, for extraction of mode shapes from spatial data. 


5 



Forsythe polynomials are defined as (Forsythe, 1957) 
P oO) = l 

Pi (x) = xp 0 (x) - a x p 0 (x) 

Pm ( x ) = x Pi ( x ) ~ cc i+1 Pi (x) - (x) 

with 

a i + 1 = Z x u fa )] 2 ; Z ( x p )J 2 

M= 1 M ' 1 

N JV 

A = Z X // -Pi (x M )Pi-l( X fi) / H [Pi ( X M )] 2 

fX=l /J=l 


( 2 . 1 ) 


( 2 . 2 ) 

(2.3) 


2.2 Estimation of adjacent mode shapes 

For an Euler beam excited at a point x = a along its length, by a force F 0 e ia * , the 
steady state response can be written, in standard form, as 

= -® s ) '((<»? -® 2 ) 2 + 4 f > ! } 

i=x (2-4) 

-i2F Wj {x)y ^co j( dI{<,co] - co 2 ) 2 +4<gV}> ,<ur 

In the above 

F = force /beam mass = F 0 / pA l 
C 0 j = yth natural frequency 
y j = yth mode shape 

= yth equivalent modal damping ratio 

Employing the Forsythe polynomials, the j th mode shape of the beam is expressed as, 

¥j( x ) = Yj S jkPk( x ) ( 2 - 5 ) 

Jfc=0 

where 

S, t j= 1,2,3...; k- 1,2,3.. J, 

are the constants to be estimated from the measured data. 
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Defining 


( 2 . 6 ) 

(2.7) 


A j = Fy/jiaXca) -a 2 ) !{{co) -co 2 ) 2 + 4% 2 co 2 } 
Bj - 2Fy/ j (a)ca jCO^j I {(a>] - co 2 ) 2 +4 %)&)} 


The displacement response of the beam, equation (2.4), can be rewritten as 

oo 

X*) = X [Aj y } (X) - iB J W] (x)]e ,a * (2.8) 

7=1 


Substitution of equation (2.8) in (2.5) gives 


x*)=S 

7=1 


(A, W) 


Jt = l 


(2.9) 


Considering that only two natural modes (with natural frequencies and co m+l ) of 
the beam, lying on either side of the excitation frequency n>(i.e. co m < co < ® m+1 ), 

predominantly contribute to the response (ignoring the contribution of the natural 
modes lying farther away from the excitation frequency), the beam response of 
equation (2.9), reduces to the approximation 

X*) = V.A m + -fS„ l )t{S„ ufl Wr (2.10) 

4=1 4=1 

The magnitudes of the real and imaginary components of the response then, are, 

nW = '>.'t^p t (x) + A., l f d S MJ ,p l ( x ) 

*« 1 1 

I'M = + (2.11) 

4=1 4=1 


Let (x^ ) and v i (x M ) be the real and imaginary components of the response 
amplitude actually measured, experimentally, at stations x M , (p. = 1, N) along the 
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length of the beam. The corresponding Forsythe form of approximation formulated 
through equation (2.10), can be written as 


h h 

Y r Op ) = X S mk p k (X M ) + A m ^Yj S m+Uk p k Op ) 

4=1 4=1 

h h 

Y , Op ) = B m Y, S mk p k (X M ) + £ m+1 £ s m+lk p k (x ) 

4=1 4=1 


The least square error between the measured amplitudes and 
approximations of equation (2.12) are 


1 

£•, =EKOp)-v' J o^)] 2 

^=i 


Noting the orthogonality property of the Forsythe polynomials 

N 

X>,Op)/>,Op)= 0 * 

P=1 


and minimizing the errors with respect to the constants S mk and S m+l k 
the following set of equations are obtained 

^m^mk ^m+l ^ m+\,k ^k 

Bm^mk + B m +\Sm+l,k ~^k \,2, h 

with 


C 4=Z V r0p)P4 0p)/Z^*0p)] 2 


N 


N 


^4=E ^0p)^0p)/Zb40p)] 2 

p=i p=i 


( 2 . 12 ) 

their Forsythe 

(2.13) 

(2.14) 

(k = l,2,....h) 

(2.15) 

(2.16) 
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The expressions for constants S mk and S m+lk can be written, from equations (2.14) 
and (2.15) as, 

S m k = ( c k B m+ 1 / A m+I ~ d k ) / [(B m+ 1 ! A m+] - B m / -4 m )T m ] 

=i c kBJA m -d k )/[(B m+l /A m+l -B m /A m )(-A m+l )] (2.17) 

The orthogonality properties of the mode shapes of a beam of constant properties over 
its span (x=0 to x= /) require that 

i 

jV,( x )^( x ) = ° * (2. 1 8) 

0 

l 

JV< (x) 1 /' j(x) = 0 i*/ (2.19) 

0 

In terms of the Forsythe approximation the above modal orthogonality implies, 

(V h h 

ZE S*P k (x M )Z 

M=\ k = 0 Jfe=0 

iV A h 

Z 'L s *p" k ( X „)Z * ( x „) = ° * (2-20) 

//=l £=0 £=0 

Substitution in the above, from equations (2.17) for modes m and ttc+ 1 and noting the 
relations (2.18) and (2.19), one obtains 

ZZ = 0 (2.21) 

//*l A:=0 

Z Z ( C * B rn + 1 / )A ( X „ )Z ( C * B m f A m~ d k )Pk ( X „ ) = 0 (2-22) 

p=l i=0 *=0 

Using the values of c k ,d k (computed from experimental data, equation (2.16)) 
equations (2.21) and (2.22) can be readily solved for the ratios B m / A m and 
B m+ , / A m+l . These ratios are further substituted in equations (2.17) to obtain the 
constants and S m+X k as multiples of ^ m and A m+1 respectively. 
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Usage of S mk and S m+lk , thus determined, in equation (2.17), give the 
approximations for the mth and (m+l)th mode shapes, y/ m (x) and ¥m+ i( x )> 
respectively. 

2.3 Natural Frequency Estimates 

The mode shapes of a beam, irrespective of its boundary conditions, are governed by 
the differential equation, 

'f'""«(x) = X 4 mV'm ( x ) (2.23) 

where ( ) denotes differentiation with respect to x. 

The spatial frequency is related to the natural frequency co m of the beam as, 
xl=cDl(pA/EI) (2.24) 

{p, A, E, / being the mass density, area of cross section, modulus of elasticity and area 
moment of inertia of beam). 

The spatial frequency can be obtained from (2.23) as, 

zt = S mk pi (X M )/S mjkPk (X M ) (2.25) 

and the natural frequency co m can be readily calculated from equation (2.24). 

2.4 Remarks 

Participation of two modes, lying on either side of the excitation frequency, in the 
overall response of the beam, can be approximately quantified through the procedure 
developed. Successive application of the algorithm at different excitation frequencies 
for estimation of desired number of mode shapes and their refinement is described in 
the next chapter. 
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CHAPTER 3 


SUCCESSIVE MODAL APPROXIMATION 


The procedure developed for estimation of adjacent modal participation, in the 
previous chapter, is employed to successively estimate the characteristics of the 
desired number of modes. The step-by-step procedure to carry out these 
approximations is described here, along with numerical illustrations for some of the 
common type of boundary configurations, namely, simply-supported, clamped- 
clamped and cantilever. The accuracy of the estimates and mode shape refinement is 
also discussed. 

3.1 Algorithm for successive modal approximation 

Following the procedure described in the previous chapter, the characteristics of 
desired number of modes can be obtained by repeated application of the algorithm, 
through appropriate selection of 

(i) frequency of excitation force and 

(ii) location of the excitation force. 

In a practical situation the boundary conditions of beams, may not behave ideally, i.e. 
the boundary configuration may not ideally be simply supported or fixed-fixed or 
cantilever etc. Since the focus of the present study is to obtain an accurate mode shape 
description, one of its major application can be in the area of obtaining correct 
definitions of the boundary conditions of a^given beam configuration. Though, this 
task has not been attempted in the present study, it has been kept as a guideline in 
developing the step-by-step procedure for successive modal estimation. 

Depending on boundary configurations a beam can be expected to have natural modes 
from any one of the following categories- 

(i) nodes at both ends (e.g. simply-supported, fixed-fixed, one end fixed- - 
other end simply-supported etc.) 
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(ii) node at one end and anti-node at the other end (e.g. cantilever i.e. fixed at one 
end and free at the other end) 

(iii) anti-nodes at both the ends (e.g. free-free case) 

(Beams with intermediate supports are not discussed here) 

The first mode in category (iii) corresponds to the zero natural frequency case, which 
signifies rigid body mode. The higher modes, however, belong to category (i), and 
therefore the vibration mode categories reduce to the first two listed above. Modal 
approximation algorithms are developed for Category (i) and Category (ii) separately. 

Category (i) nodes at both ends (e.g. simply-supported, fixed-fixed, etc.) 

Typical mode shapes of simply supported and fixed-fixed beam configurations are 
shown in Figures 3.1 (a) and (b), respectively. The mode shapes can be seen to be of 
two types - 

(a) symmetric with respect to the mid-point (odd modes, 2/ - 1 , i = 1,2 ) 

(b) anti-symmetric with respect to the mid-point (even modes, 2 i i = 1,2 ) 

Step 1 : Estimation of Symmetric Modes 

The first step is to select the location of the excitation force at, approximately, the 
mid-point of the beam. Such a mid-point excitation is expected to primarily energize 
the symmetric modes (the nodal points of anti-symmetric modes being at the beam 
mid-point). The frequency of the excitation force needs to be chosen appropriately. 
An a-priori estimate of the fundamental frequency, if available, from basic test tests, 
like the rap-test, can help (the case, of no a-priori knowledge of the fundamental 
frequency, is discussed later). In case such an estimate, p l , is available, the excitation 
frequency, <y, , can be chosen to be slightly above this frequency (say co x = \2p x ) so 
that the response is primarily contributed from the first two odd modes (i.e. mode 1 
and mode3). Ignoring the contributions from higher symmetric modes (mode 5, mode 
7 etc.), the approximations for the first and third mode shapes can be obtained, 
through the methodology discussed in the previous chapter. Simultaneously, 
approximations of the first and third natural frequencies are also obtained. 
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Fig. 3.1(a) Exact mode shapes of simply supported beam 
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Fig. 3.1(b) Exact mode shapes of fixed-fixed beam 





Armed with an approximation of the third natural frequency, p 3 , the next step would 
be to excite the beam, at its mid-point, at a frequency, co 2 , slightly above its third 
natural frequency (say a> 2 =l.2p 3 ). The response now, has contributions from 

primarily from the third and fifth modes. However, it may also contain significant 
contribution from the first mode. Since, the first mode has already been approximated, 
its contribution can be subtracted from the response and the resultant can be fed as 
input to the algorithm for estimation of the third and fifth mode shapes and natural 
frequencies. The next excitation frequency, a> 3 , can be above the fifth natural 

frequency p s and the process can be continued further for estimation of the desired 

number of symmetric modes. All these estimates can be termed as ‘first 
approximations’. 

Step 2: Estimation of Anti-Svmmetric Modes 

The ‘first approximations’ of the anti-symmetric modes can be obtained through a 
similar step-by-step procedure. From the knowledge of the nodal points of the 
symmetric modes, estimated as described above, the location of excitation force, can 
be chosen close to such a nodal point so as to suppress the contribution of the 
corresponding mode in the response. Choosing the nodal point of mode 3 (at a 
distance 1/3 or 2//3, from one of the ends), the beam can be given an excitation at a 
frequency, co 4 , such that p ] < co 4 < p 3 . The modes that will get excited are 1, 2, 4, 

5,... The contributions of mode 1 and 5, which have already been estimated above, 
can be subtracted from the total response and the resultant can be fed as input to the 
algorithm for obtaining the ‘first estimates’ of the anti-symmetric modes 2 and 4. The 
process can be continued further, through excitation at the nodal points of the 
symmetric mode 5 and so on, to get the ‘first estimates’ of the higher anti-symmetric 
modes. 

Step 3: Refinement 

It should be noted that, the algorithm gives an estimate of two modes, at a time. The 
‘first estimates’ carry errors, since the algorithm treats the beam response to be 
constituted of just the two modes that it estimates. Refinement can be carried out, by 
repeating Steps 1 and 2. Since the ‘first estimates’ of the modes are now available, the 
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contribution of all others, except for the two modes being refined can be subtracted 
from the total response of the beam, before feeding it as input to the algorithm. The 
iterations can continue in this manner, till the difference between two successive sets 
of estimates is less than desired accuracy. 

Category (ii) node at one end - anti-node at the other end (e.g. cantilever) 

Typical mode shapes of a cantilever beam are shown in Figures 3.1 (c). Unlike, the 
previous category, the modes of a cantilever beam cannot be classified as symmetric 
or anti-symmetric and it would be simpler to keep the location of excitation 
identically at the free-end, for all runs of the algori thm 

First Estimates: 

As in the previous category, if an a-priori knowledge of the first natural frequency, 
Pi , exists the excitation frequency, co x , can be chosen to be slightly above p x (say, 
co x = 1.2/?, ), so that the primary contributions to the beam response are from mode 1 
and 2. ‘First estimates’ of these two modes can be obtained by feeding this response 
as input to the algorithm. Having estimated the natural frequencies of modes 1 and 2, 
the excitation frequency can be slightly above the second natural frequency (say, 
co 2 =\.2p 2 ). The contribution of mode 1, through its ‘first estimate’, can be 
subtracted from the beam response, so that the resultant has contributions mainly, 
from modes 2 and 3. Feeding this resultant as input to the algorithm, modes 2 and 3 
are determined. The procedure can be continued by, next, exciting the' beam at a 
frequency slightly higher than the highest natural frequency determined, for obtaining 
the ‘first estimates’ of the higher modes. 

Refinement: 

The ‘first estimates’ can be refined, by repeating the above process. However, this 
time the contributions of all modes, except the two of interest, are to be subtracted 
from the beam responses, through their ‘first estimates’. The resultant of the response, 
thus, is a better representation of the contributions of the two modes being estimated 
by the algorithm. 
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Fig. 3.1(c) Exact mode shapes of cantilever beam 



In the procedures described above, the starting point was excitation of the beams at a 
frequency slightly above their fundamental frequency. In cases where an a-priori 
knowledge of this fundamental frequency may not exist, the frequency chosen to 
excite the beam lies between the z'th and (z'+l)th natural frequencies, the algorithm 
takes the entire response of the beam to be constituted of the z'th and (z'+l)th modes 
and gives estimates of these two mode shapes. The numerical values of z and (z'+l) can 
be readily found out from the mode shapes, by counting the number of nodes and anti- 
nodes. The algorithm will also give estimates of the z'th and (z+l)th natural 
frequencies. By judiciously choosing the next excitation frequency lower than the fth 
natural frequency and continuing the process, the fundamental frequency can be 
traced back. 

It should be noted here, that the procedures outlined can not be employed as a black 
box. Step-by-step analysis of the results is essential along with the step-by-step 
application of the algorithm. 

The procedures have been illustrated through computer simulation. The illustrations 
are carried out in terms of nondimensional parameters. 

3.2 Non-Dimensional Parameters 

The response of the beam, subjected to a desired excitation, is simulated through a 
computer program employing standard response formulae (Timoshenko, 1928). This 
simulated response is treated as the ‘measured’ response of the beam and fed as input 
to the modal estimation algorithm described in the previous chapter. 

The response of a beam, subjected to a force F 0 e“* at a point x= a, along its length, 
can be written in nondimensional terms as (refer, equation 2.4) 

= V,(X) V“ (3-1) 

where 

H (x) = y(x) / y st 
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and 


y s t = F I pAlco\ 

r } - co / go t 

Aj = y j (flXl - r; ) / ((1 - r; ) 2 + 4 t;rj ) 

B, = Vl (d)(2r^l((\-r]f +4 $rj) (3.2) 

The geometric and material properties of the beam are incorporated in a 
nondimensional parameter, k , defined as 

k = [EI / pgA! 3 f /2) (3.3) 

( g , the gravitational acceleration is taken as 9.81 m/s 2 ). 

The nondimensional parameter, k , establishes the relationship between the yth spatial 
frequency Xj and the yth natural frequency co } , through the following expression 

(refer equation 2.24) 

0> m = Z 2 m *(glY 2 ( 3 - 4 ) 

Referring to equations (2.4)-(2.9) and taking the response to be constituted primarily 
of the mth and (/n+l)th beam modes, which lie on either side of the excitation 
frequency, co , the magnitudes of the real and imaginary components of the response 
the nondimensionalised Forsythe approximation of the response are 

h h 

H r (x) = A m 'Z S mkPk (x) + A m+1 2 S m+lk p k (. x ) 

1 *=1 

h _ h 

H t (x) = Bm YS mkPk (*) + B m+l J'jS /n+l,jfc Pk (*) (3-5) 


Defining, the real and imaginary components of the measured response in 
nondimensional terms as 
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q r (*) = 1 ( F / P Al ) 

q l {x)=v i {x)co\ / (F / pAl) 


(3.6) 


the non-dimensionalised least square errors are 

l =£[#,« -9,W] 2 
//=1 

<=Xww-«,wr o-v 

The non-dimensionalised errors are now minimiz ed with respect to 
S mk and S m+l k (k = 1,2 h ) for estimation of the constants and S m+[ k and 

subsequently the mth and (m+l)th modal parameters, as described in the previous 
chapter. 

3.3 Numerical illustrations 

The numerical studies for successive modal parameter estimation has been carried out 
for the following three beam configurations. 

(i) Simply-Supported 

(ii) Cantilever (Clamped-Free) 

(iii) Fixed-fixed (Clamped-Clamped). 

The estimation, using the computer simulated response is carried out for a value of the 
nondimensional beam parameter, k= 0.1. The procedure is illustrated for estimation 
of the first five beam modes. The damping ratios, , are chosen as 0.02 for 

illustration. The nondimensional input force amplitude, Ft pgAl is kept as 100 in all 
the cases. 

Simply supported: 

The illustration is carried, with an assumption that an a-priori knowledge of the 
fundamental frequency exists. (This frequency can be readily obtained, as mentioned 
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earlier, through simple rap-tests and analysis of the consequent transient response) 

The fundamental frequency, in this case (with k = 0.1) is 3.09 rads/s, (Rao-1986) 

The step-by-step modal estimation procedure is implemented as follows. 

(a) The beam is harmonically excited at its mid-point at a frequency of 3.71 rads/s 
(= 1.2 x 3.09 rads/s). Figures 3.2 (a) and (b) show the in-phase and out-of- 
phase components of the computer simulated spatial response H(x), using 
equation (3.2). The response is plotted as a function of the nondimensional 
distance, x/l, along the beam. The response of the beam is constituted of 

contributions from all the symmetric natural modes (modes 1, 3, 5, ) of 

the beam. The estimation procedure is applied by picking up this spatial 
response at hundred number of points along the beam. Reference can be made 
to a previous study (Vaze, 1998) for a discussion on he influence of the 
number of these measurement points and the influence of possible random 
measurement errors. The algorithm was found to be robust in the presence of 
measurement noise. The estimates of modes lying adjacent to the excitation 
frequency (mode 1 and mode 3, in this case) are shown in Figs. 3.3 (a) and (b), 
along with the exact mode shapes of the beam. The estimate of the 
corresponding spatial frequencies, % , are shown in Figs. 3.4 (a) and (b) and 
the error in mode shapes is as shown in Fig. 3.5 (a) and (b). It should be noted 
here that estimates of the spatial frequencies, x ■> (Rao-1986) can be obtained 

through application of equation (2.25) at any point x M , along the beam. Figs. 
3.4 (a), (b) show the estimates of x at a11 measurement points. The exact 
spatial frequencies are also shown in the graphs as straight lines. The error in 
the estimates of x can be seen to be considerable. Also, the error is not 
identical across the beam length. It needs to be emphasized that the 
computation of the spatial frequency % » involves a fourth order numerical 
differentiation. A central-difference scheme is employed for this numerical 
differentiation, which carries inherent errors. 
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Fig 3.2(a) In-phase response of simply supported beam 
(Location of excitation force, xjl- 0.5 
Excitation frequency = 3.71 rad/s) 
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Fig. 3.2(b) Out of phase response of simply supported beam 
(Location of excitation force, xjl= 0.5 
Fynitatinn fr^mipnnv = 3 71 raHM 
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Fig. 3.3(a) First approximation of mode 1: simply supported beam 
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Fig. 3.3(b) First approximation of mode 3: simply supported beam 
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It was observed, during the present study, that the numerical differentiation 
scheme is not ideally suited for the spatial frequency and consequent natural 
frequency estimation. It can be said that natural frequency estimation should 
be carried out through conventional modal analysis routes (free or forced 
vibration tests), while the focus of spatial analysis of the present study should 
remain confined to obtaining accurate descriptions of mode shapes. However, 
for the present work, the means of the % curves are taken as the estimates of 
the spatial frequencies. The corresponding natural frequencies, consequently 
computed using equation (2.24) are 3.081 rads/s. for mode 1 and 32.480 
rads/s. for mode 3. 

(b) The excitation is next applied at the mid-point, at a frequency of 38.856 rads/s. 
(=1.2 x 32.380 rads/s.). The beam response is shown in Figs. 3.6 (a), (b). The 
contribution of the first mode, estimated above, is subtracted from this 
response. The resultant, shown in Figs. 3.7 (a), (b), is fed as input to the 
algorithm and modes 3 and 5 are estimated. The estimated mode shapes, 
spatial frequencies and the error in mode shapes are given in Figs. 3.8 (a), (b), 
3.9 (a), (b) and 3.10 (a), (b) respectively. The natural frequencies computed 
using equation (2.24) are 31.528 rads/s. for mode 3 and 37.409 rads/s. for 
mode 5. 

(c) The anti-symmetric modes are estimated now, by choosing the excitation force 
location at x/l =1/3, which is a nodal point of mode 3. The simulated response 
is given in Figs. 3.11 (a), (b). The excitation frequency is the mean of the 
estimated first and third natural frequencies estimated above (17.781 rads/s. 
= {3.08 1 + 32.480 }/2 ). The beam response is shown in Figs. 3.1 1 (a), (b). The 
contribution of the first mode, estimated in (a), is subtracted from this 
response. The resultant, shown in Figs. 3.12 (a), (b), is fed as input to the 
algorithm and modes 2 and 4 are estimated. The estimated mode shapes, 
spatial frequencies and the errors in mode shapes are given in Figs. 3.13 (a), 
(b), 3.14 (a), (b), and 3.15 (a), (b) respectively. The natural frequencies 
computed using equation (2.24) are 16.905 rads/s. for mode 2 and 75.386 

rads/s. for mode 4. 
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Fig. 3.7(b) Out of phase response of simply supported beam after subtractin 

the contribution of 1 st mode 
(Location of excitation force, x/l= 0.5 
Excitation frequency = 38.856 rad/s) 
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Fig. 3.7(a) In-phase response of simply supported beam after subtracting 
the contribution of 1 st mode 
(Location of excitation force, x//= 0.5 
Excitation frequency = 38.856 rad/s) 
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Fig. 3.9(a) First approximation of Spatial frequency: mode 3: simply supported beam 
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.9(b) First approximation of Spatial frequency: mode 5: simply supported beam 
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Fig. 3. 1 0(a) Error in the first approximation of mode 3; simply supported beam 
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Fig. 3 . 1 1 (a) In-phase response of the simply supported beam 
(Location of excitation force, x/l= 0.33 
Excitation frequency = 17.781 rad/s) 
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Fie. 3.1 Kb) Out of phase response of the simply supported beam 
(Location of excitation force, x/l - 0.33 





Displacement 


0.4 
0.2 

-t— * 

c 
© 

E 
© 
o 
J5 

Q. 
w 

b 

- 0.6 
- 0.8 
-1 

Length along the beam 

Fi = . 3.12(a) In-phase response of simply supported beam after subtracting 

the contribution of 1 51 mode 
(Location of excitation force, xjl = 0.33 
Excitation frequency = 17.781 rad/s) 

2.5 
2 

1.5 
1 

0.5 
0 

-0.5 
-1 
-1.5 
-2 
-2.5 

0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 

Length along the beam 

Fig. o. 12(a) Out of phase response of simply supported beam after subtracting 

the contribution of 1 st mode 
(Location of excitation force, xjj- 0.33 
Excitation frequency - 17.781 rad/s) 

32 











Spatial frequency ' x 2 ’ (db) Spatial frequency ‘ 1 (db) 



Length along the beam 

Fig. 3.14(a) First approximation of Spatial frequency: mode 2: simply supported beam 



Length along the beam 


Fig. 3. 14(b) First approximation of Spatial frequency: mode 4: simply supported beam 





The above forms ‘the first estimates’ of the modal characteristics of the simply 
supported beam under consideration. Refined mode shapes and spatial frequencies are 
obtained by carrying out repeated iterations of the procedure, in a manner described 
previously. The refinement obtained in the first five mode shapes are shown in Figs. 
3.16 (a)-(e), while the corresponding data on the spatial frequencies is depicted in 

Figs 3.17 (a)-(e). The natural frequencies co m , obtained in each iteration, through the 

formula given in equation (2.24), are listed in Table 3.1, along with the corresponding 
exact values. 


Table 3.1 Estimated and Exact Natural Frequencies 
of Simply Supported Beam 



Mode 1 

Mode 2 

Mode 3 

Mode 4 

Mode 5 

Iteration 1 

wmmm 

16.905 

32.480 

75.386 


Iteration 2 

i 

12.638 

30.597 

54.848 

IH339H 

Iteration 3 


12.607 

29.572 

51.298 

72.651 

Exact 

3.091 

12.365 

27.821 


77.281 


Fixed -Fixed ( clamped-clamped): 

The procedure in this case, follows a pattern similar to that in the previous case, since 
both configurations carry antinodes at the ends and belong to the same category 
defined earlier. The fundamental frequency, in this case, with non-dimensional 
parameter, k = 0.1 was found out to be 7.0rads/s. 

(a) The beam is harmonically excited at its mid-point (x/l= 0.5), at a frequency 
equal to 8.4 rads/s (= 1.2 x 7.00rads/s), so that the response is primarily 
constituted of the fundamental mode and the next symmetric mode (i.e. mode 
3). Figures 3.18 (a) and (b) show the ‘first estimates’ of modes 1 and 3. The 
‘first estimates’ of the corresponding spatial frequencies, % , are shown in 
Figs. 3.19(a), (b) along with the exact values. The natural frequencies 
computed from the mean of the % curves, are 7.151 rads/s. for mode 1 and 


.42.226 rads/s. for mode 3. 
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Fig. 3.16(c) Refinement in mode shapes: mode 3: simply supported beam 
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Fis. 3.16(d) Refinement in mode shapes: mode 4: simply supported beam 
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Fig. 3.17(a) Refinement in Spatial frequency: mode 1: simply supported beam 
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Fig. 3.17(b) Refinement in Spatial frequency: mode 2: simply supported beam 
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Fig. 3.17(d) Refinement in Spatial frequency: mode 4: simply supported beam 
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Fig. 3.19(a) First approximation of Spatial frequency: mode 1: fixed-fixed beam 
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Fig. 3.19(b) First approximation of Spatial frequency: mode 3: fixed-fixed beam 




(b) For estimation of the next set of symmetric modes (3 and 5), mid-point 
excitation is given, to the beam, at a frequency of 50.671 rads/s. (=1.2 x 
42.226). The contribution of the first mode, estimated in step (a), is subtracted 
from this response. The resultant is fed as input to the algorithm and modes 3 
and 5 are estimated. The estimated mode shapes and spatial are given in Fig. 
3.20 (a), (b), Fig.3.21 (a), (b) respectively. The natural frequencies computed 
from the mean of the spatial frequency curves are 33.1 10 rads/s. for mode 3 
and 72.951 rads/s. for mode 5. 

(c) In order to estimate the anti-symmetric modes (2,4, ), excitation is provided 

at the nodal point (x/7=0.35) of the third mode estimated in step (b). The 
excitation frequency is 24.688 rads/s. (=mean of 1 st and 3 rd estimated natural 
frequencies). The contribution of the first mode, estimated in (a), is subtracted 
from the response. The resultant is fed as input to the algorithm and modes 2 
and 4 are estimated. The estimated mode shapes and spatial frequencies are 
given in Figs. 3.22 (a)-(b) and 3.23(a)-(b), respectively. The estimates for the 
natural frequencies computed using the mean of the spatial frequency curves 
are 16.693 rads/s. for mode 2 and 41.097 rads/s. for mode 5 respectively. 

Refinements of the above ‘first estimates’ are carried out in a manner similar to that 
described in the simply supported case. This refinement in mode shapes and spatial 
frequencies is illustrated in Figs. 3.24 (a)-(e) and Figs 3.25 (a)-(e), respectively. The 

natural frequencies a > m , obtained through the formula given in equation (2.24) are 
listed in Table 3.2, for each iteration, along with the corresponding exact values. 


Table 3.2 Estimated and Exact Natural Frequencies 
of fixed-fixed beam 



Model 

Mode 2 

Mode 3 

Mode 4 

Mode 5 

Iteration 1 

7.150 

16.693 

42.226 

41.097 

72.951 

Iteration 2 

6.799 

17.908 

41.094 

52.221 

84.067 

Iteration 3 

6.799 

17.978 

40.643 

52.671 

84.094 

Exact 

7.000 

19.316 

37.868 

62.597 

93.510 1 
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Fig. 3.22(a) First approximation of mode 2 : fixed-fixed beam 
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Fig. 3.22(b) First approximation of mode 4: fixed-fixed beam 
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Fig. 3.25(a) Refinement in Spatial frequency: mode 1 : fixed-fixed beam 
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Fig. 3.25(b) Refinement in Spatial frequency: mode 2: fixed-fixed beam 
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Fig. 3.25(e) Refinement in Spatial frequency: mode 5: fixed-fixed beam 



Cantilever fClamped-free - ): 

The fundamental frequency, for the cantilever (for /c = 0.1) is found out to be 
l.lOrads/s. 

(a) First the beam is harmonically at its free end at a frequency of 1.32 rads/s 
(=1.2 x 1.1 rads/s). The computer simulated spatial response H(x), using 
equation (3.2), is given in Figs. 3.26 (a) and (b). The response in this case is 

constituted of contributions from all the natural modes (modes 1, 2, 3, 4 ) 

of the beam. For the first iteration, making an assumption that the modes (1 
and 2) lying closest to the excitation frequency, are the major contributors to 
the overall response, the algorithm is implemented to obtain approximations 
for the mode shapes (Figs. 3.27 (a)-(b)) and spatial frequencies (Figs. 3.28 (a)- 

(b)) of modes 1 and 2. The corresponding natural frequencies obtained from 
the means of the x curves, are 0.973 rads/s. for mode 1 and 6.221 rads/s. for 
mode 2. 

(b) Approximations are obtained next for the set constituted of mode numbers 2 
and 3, by choosing the excitation frequency above the second natural 
frequency estimated above 7.465 rads/s. (=1.2 x 6.221 rads/s). The 
contribution of the first mode, estimated in step (a), is subtracted from the 
response and fed as input to the algorithm. The estimated mode shapes and 
spatial frequencies are shown in Figs. 3.29(a)-(b), 3.30(a)-(b), respectively. 
The natural frequencies obtained are 5.785 rads/s. for mode 2 and 21.365 
rads/s. for mode 3. 

(c) The process is continued in a similar manner to obtain approximations for 
successive sets of modes. The excitation is always given at the free end at a 
frequency higher than the highest estimated in a previous run. Approximations 
have been obtained up to mode 5. 

The ‘first estimates’ are refined, by iterating through the above procedure. Since the 
‘first estimates’ of all the modes of interest are available, the possible participation of 
all modes not lying adjacent to the selected excitation frequency can be computed and 
subtracted from the overall response. The resultant thereby forms the refined input 
constituted of the participation of only those modes, which are being estimated by the 
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Fig. 3.26(a) In-phase response of the cantilever beam 
(Location of excitation force, xjl~ 1.0 
Frequency of excitation = 1.32 rad/s) 
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Fig. 3.26(b) Out of phase response of the cantilever beam 
(Location of excitation force, x\l = 1 .0 
Frequency of excitation = 1.32 rad/s) 
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Fig. 3.27(a) First approximation of mode 1; cantilever beam 



Fig. 3.27(b) First approximation of mode 2; cantilever beam 
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Fig. 3.28(a) First approximation of Spatial frequency: mode 1: cantilever beam 















algorithm. The input to be fed to the algorithm is refined by computed by subtracting 
the contributions of the modes not of interest from the overall beam response. The 
refinement obtained in the first five mode shapes is shown in Figs. 3.31(a)-(e), while 
the corresponding data on the spatial frequencies is depicted in Figs 3.32 (a)-(e). The 

natural frequencies co m , obtained through the formula given in equation (2.24) are 
listed in Table 3.3, for each iteration, along with the corresponding exact values. 


Table 3.3 Estimated and Exact Natural Frequencies 
of cantilever beam 



Mode 1 

Mode 2 

Mode 3 

Mode 4 

Mode 5 

Iteration 1 

eesoh 

6.221 

21.365 

42.102 


Iteration 2 

K! ■ 

6.823 

18.322 

41.030 


Exact 

1.101 

6.955 

19.320 

37.867 



3.4 Remarks 


The results obtained through the numerical illustration are self-explanatory. The 
procedure can be seen to obtain very accurate estimates of mode shapes. For example, 
the maximum error in the simply supported case ranges from about 0.015 % for mode 
1 to 10 % for mode 5. A similar trend can be observed in the other beam 
configurations. Presently, the estimation is restricted to the fifth mode in all cases. For 
improvement in the fifth mode approximations, higher modes (6 th , 7 th ,...) need to be 
considered in the analysis. The errors are found to be maximum in the vicinity of the 
nodal points. The exact displacement at these points is zero. The computer algorithm, 
however do not approximate these displacements as exact zeros. The same trend is 
observed in the approximations of the spatial frequencies. The errors at the nodal 
points is more pronounced in the case of spatial frequencies since a fourth order 
numerical derivative is involved in the computational procedure. 
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Fig. 3.31(c) Refinement in mode shapes: mode 3: cantilever beam 
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Fig. 3.31(e) Refinement in rr 









Spatial frequency 4 ’ (db) Spatial frequency 4 % ’ ( db ) 



0.1 0.2 0.3 

Length along the beam 

Fig. 3.32(c) Refinement in Spatial frequency: mode 3: cantilever beam 



Fig. 3.32(d) Refinement in Spatial frequency: mode 4: cantilever beam 








CHAPTER 4 


CONCLUSIONS 


The work presented in this thesis is an extension of the existing Spatial Modal 
Analysis procedures, dealing with estimation of single (or two) natural modes, of a 
vibrating beam, which is (are) close to the frequency of excitation. An algorithm has 
been developed, during the present study, for step-by-step approximation of 
successive beam modes. The procedure involves variation of the frequency and 
location of the excitation force. The general procedure constitutes of estimation of 
mode shapes, as well as calculation of natural frequencies, from steady state spatial 
data of beam excited at a specified frequency and location. The deflected beam shape 
is approximated by employing recursive form of Forsythe polynomials and Least 
Square technique is used to curve fit the measured responses and obtain the 
coefficients of Forsythe polynomials. Orthogonality properties of the normal mode 
functions are further employed to decompose the Forsythe polynomial representation 
into the individual mode shapes of the beam contributing to the response. The 
individual mode shape estimates are further processed to fit the governing fourth 
order spatial equation of beams to obtain the spatial modal frequency. Successive 
beam modes are approximated, step-by-step, by changing the frequency and location 
of excitation force. The procedure has been illustrated through computer simulation 
for beams with standard boundary configurations. The procedure is found to give 
accurate descriptions of the mode shapes. The estimates of the natural frequencies 
however, are not accurate. The natural frequencies have been computed through the 
estimates of the spatial natural frequencies, involving fourth order numerical 
differentiation. The central-difference scheme employed for numerical differentiation, 
carries inherent errors. It can be said that the present spatial modal analysis scheme is 
not ideally suited for natural frequency estimation, which is carried out more 
accurately through conventional modal analysis routes (free or forced vibration tests). 
The focus of the spatial analysis procedures, of the present study, should remain 
confined to obtaining accurate descriptions of mode shapes. 



The procedure needs to be experimentally validated. Since spatial modal analysis 
procedures require intensive response data from the space coordinates, conventional 
contact type sensors like accelerometers, which can provide data at limited number of 
spatial locations, are not suitable for such applications. A non-contact type of sensor 
like the Laser Vibrometer is essential to get the dynamic information at sufficient 
number of stations on the beam. Influence of the number of Forsythe polynomials in 
the mode shape approximations, also needs to be investigated. An investigation is also 
required to explore and account for the errors incurred due to inaccuracies in the force 
application point. The study also opens the possibility of developing similar 
procedures for plate vibration applications. Possibilities of application of multiple- 
point excitations can also be investigated. 
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